How to Measure Subdiffusion Parameters 
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We propose a method to measure the subdiffusion parameter a and subdiffusion coefficient D a 
which are defined by means of the relation {x 2 } = r fi+ct) ^ wriere {x 2 ) denotes a mean square 
displacement of a random walker starting from x — at the initial time t — 0. The method exploits 
a membrane system where a substance of interest is transported in a solvent from one vessel to 
another across a thin membrane which plays here only an auxiliary role. We experimentally study 
a diffusion of glucose and sucrose in a gel solvent, and we precisely determine the parameters a and 
D a , using a fully analytic solution of the fractional subdiffusion equation. 



PACS numbers: 05.40.-a, 66.10.-x 

Subdiffusion occurs in various systems. We mention 
here a diffusion in porous media or charge carriers trans- 
port in amorphous semiconductors [1,2]. The subdiffu- 
sion is characterized by a time dependence of the mean 
square displacement of a Brownian particle. When the 
particle starts form x — at the initial time t = this 
dependence in a one-dimension system is 



<* 2 > 



2D a 



r(i 
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where D a is the subdiffusion coefficient measured in the 
units [m 2 /s a ] and a obeys < a < 1. For a = 1 one 
deals with the normal or Gaussian diffusion character- 
ized by the linear growth of {x 2 } with t which results from 
the Central Limit Theorem applied to many independent 
jumps of a random walker. The anomalous diffusion oc- 
curs when the theorem fails to describe the system be- 
cause the distributions of summed random variables are 
too broad or the variables are correlated to each other. 
The subdiffusion is related to infinitely long average time 
that a random walker waits to make a finite jump. Then, 
its average displacement squared, which is observed in a 
finite time interval, is suppressed. 

The subdiffusion has been recently extensively stud- 
ied, see e.g. [1-4]. While the phenomenon is theoreti- 
cally rather well understood there are very a few exper- 
imental investigations. There is no effective method to 
experimentally measure a and D a . In the pioneering 
study [4], where D a w&s determined experimentally for 
the first time, the interdiffusion of heavy and light wa- 
ter in a porous medium was observed by means of NMR. 
D a was found, using the special case a = 2/3 solution of 
the subdiffusion equation. The procedure is neither very 
accurate nor of general use. 

Our aim here is to present a method to precisely mea- 
sure a and D a . The method is described in detail in 
[5], here we give a brief account of it. For practical rea- 
sons, we choose for the experimental study a membrane 
system containing two vessels with a thin membrane in 
between which separates the initially homogeneous so- 
lute of the substance of interest from the pure solvent. A 



schematic view of the system is presented in Fig. 1. The 
membrane does not affect values of investigated param- 
eters. Instead of the mean square displacement (1), our 
method refers to the temporal evolution of the thickness 
6 of the so-called near-membrane layer which is defined 
as a distance from the membrane where the substance 
concentration C(x,t) drops n times with respect to the 
membrane surface i.e. 



C{5,t) = KC(0+,t), 



(2) 



where x = is the position of a thin membrane and k 
is an arbitrary number smaller than unity; we used k — 
0.12, 0.08 and 0.05. 
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1. Schematic view of the membrane system under 



In our previous paper [6], we demonstrated that 8(t) = 
A\ft for the normal diffusion. Studying experimentally 
the diffusion of glucose and sucrose in a gel solvent, we 
show here that 6(t) = At 1 with 7 < 0.5. A gel is built of 
large and heavy molecules which form a polymer network. 
Thus, the gel water solvent resembles a porous material 
filled with water. Since a mobility of sugar molecules 
is highly limited in such a medium the subdiffusion is 
expected. 
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For each measurement, we prepared two gel samples: 
the pure gel - 1.5% water solution of agarose and the same 
gel dripped by the solute of glucose or sucrose. The con- 
centration of both sugars in the gel was fixed to be either 
0.1 [mol/dm 3 ] or 0.07 [mol/dm 3 ] but our results appear 
to be independent of the initial concentration. The two 
vessels of the membrane system were then filled with the 
samples and the (slow) processes of the sugar transport 
across the membrane started. Since the concentration 
gradient was in the vertical direction only, the diffusion 
is expected to be one-dimensional. We used an artificial 
membrane of the thickness below 0.1 mm. The mem- 
brane was needed for two reasons. It initially separated 
the homogeneous sugar solute in one vessel from the pure 
gel in another one. It also precisely fixed the geometry 
of the whole system. 
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FIG. 2. The experimentally measured thickness of the 
near-membrane layer 5 as a function of time t for glucose with 
k = 0.05 (□), k = 0.08 (o), k = 0.12 (A), and for sucrose with 
k — 0.08 (0). The solid lines represent the power function 
At 0A5 while the dotted ones correspond to the function A\ft. 

The sugar concentration was measured by means of the 
laser interfcromctric method. The laser light was split 
into two beams. The first one went through the system 
parallelly to the membrane surface while the second (ref- 
erence one) went directly to the light detecting system. 
The interferograms, which appear due to the interfer- 
ence of the two beams, are controlled by the refraction 
coefficient of the solute which is turn depends on the sub- 
stance concentration. The analysis of the interferograms 
allows one to reconstruct the time-dependent concentra- 
tion profiles of the substance transported in the system 
and to find the time evolution of the near-membrane lay- 
ers which are of our main interest here. The experimen- 
tal set-up is described in detail in [7]. It consists of the 
cuvette with membrane, the Mach-Zchnder interferome- 
ter including the He-Ne laser, TV-CCD camera, and the 
computerized data acquisition system. 



When the sugar was diffusing across the membrane 
we were recording the concentration profiles in the vessel 
which initially contained pure gel. The examples of typ- 
ical interferograms and extracted concentration profiles 
are presented in [7] . The thickness of a near- membrane 
layer 8 was calculated from the measured concentration 
profiles C(x,t) according to the definition (2), and thus 
the thickness of near-membrane layer as a function of 
time was found. 

In Fig. 2 we present 8(t) for the glucose and sucrose 
of initial concentration 0.1 [mol/dm 3 ]. The analysis of 
errors, in particular those shown in Fig. 2, is described 
in [5]. For the glucose we present 8(t) for three values of 
k = 0.12, 0.08 and 0.05 while for the sucrose k = 0.08. 
As seen, the time dependence of 8 is well described by 
the power function A t 1 with the common index 7 = 
0.45. The lines representing ~ y/i are also shown for 
comparison. It is evident that the measured index 7 is 
smaller than 0.5. There arc some deviations of our data 
from A i 45 at t < 300s but our final theoretical formulas, 
in particular the power law behavior, hold in the long 
time approximation. 

Fitting the experimental data shown in Fig. 2, we 
found the universal index 7 = 0.45±0.005 and the param- 
eter A which depends on k; for glucose A = 0.091 ±0.004 
when k = 0.05, A = 0.081 ± 0.004 when k = 0.08, 
and A = 0.071 ± 0.004 when k = 0.12; for sucrose 
A = 0.064 ± 0.003 when k = 0.08. In each case \ 2 per 
degree of freedom was smaller than 1. 

The subdiffusion is described by the equation with 
fractional derivative [1,8] 

dc(x,t) _ a 1 -" d 2 c(x,t) 
dt "at 1 - dx? ' 1 ' 

which for a < 1 corresponds to an infinitely long average 
waiting time of the random walker - the physical situation 
in a gel solvent resembling the porous medium. We solve 
Eq. (3) in the region x > with the initial condition 
C(x, 0) = C for x < and C(x, 0) = for x > 0. In 
fact, we solve Eq. (3) for the Green's function G(x, t; xo) 
satisfying the initial condition G(x, t = 0; x ) = S(x—x ), 
and then, C(x, t) is calculated using the formula 

C(x,t) = J G(x,tx o )C(x ,0)dx o . (4) 

To find G(x, t; x ) we use the relation [5] 

G(x,t;x ) = / dt' J(0+ f; x ) G rcf {x, t - t'; 0+) , (5) 
Jo 

where x > while xq < 0; J(x, t; xq) is the flux associated 
with G(x, t; xq) which for x — gives the flow across the 
membrane; G Te { (x, t; xq) is the Green's function for the 
half-space system with x > and the fully reflecting wall 
at x = 0. 
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Using Eqs. (4,5), C(x,t) can be written as 



C(x,t)= / dt'W(t')G rci (x 1 t-t':0 + ). 



(6) 



where the function W(t), which equals 



W(t) 



dxo J(0 + , t: xq) C(xq, 0) . 



depends on the initial and boundary conditions. 

Since the subdiffusion equation is of the second order 
with respect to x, it requires two boundary conditions at 
the membrane. The first one assumes the continuity of 
the flux J, given by the generalized Fick's law [9], which 
flows through the membrane i.e. J(0~,t) = J(0 + ,t). 
However, there is no obvious choice of the second bound- 
ary condition. Therefore, we assume that the missing 
condition is given by a linear combination of concentra- 
tions and flux i.e. 



b 1 C{0-,t) + b 2 C(0+,t) + b 3 J{0,t) = 0. 
Two boundary conditions 



and 



J(0,i) = X(C(Q-,t) -G(0+,i)) 



(7) 



(8) 



(9) 



discussed in [6,10] and [11], respectively, are of the gen- 
eral form (7). The parameters a and A control the mem- 
brane permeability [6,10,11]. The adopted initial condi- 
tion combined with Eq. (7) provide 



W(t) = Co 

OO 



i 



h - b 2 i 1 -"/ 2 
d k 



(10) 



^ r(a/2 - fc(l - a/2)) i*(i-«/2) 



fe=0 

where d = b 3 ^/D a / (pi — b 2 ). The Green's function G re f, 
which enters Eq. (6), can be easily obtained by means of 
the method of images [1] as G re f (x, t; xq) = Gq(x, t; xq) + 
Go(— x, t; xq) with the known Green's function Go for the 
homogeneous system [1] . Having the explicit functions W 
and G re f, we write down, using Eq. (6), the concentration 
profile as 



C(x,t) = f dt'W(t-t') 
Jo 



(11) 



x — H 

ax 
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where H denotes the Fox function. 

We first consider the long time approximation of the 
formula (11) which corresponds to the small s limit of 



the Laplace transform L{f(t)} = L°° f(t)e~ st . Taking 
into account only the leading contribution in the small s 
limit, Eq. (11) gets the form 



C(x,t) 



2G 6i 



(6i - b 2 )a 



1 1 

^ 



(12) 



The solution (12) can be also obtained directly from 
Eq. (11), taking into account only the k = term in 
the expansion (10). 

The series (10) can be approximated by the first term 
if d <C t 1- "/ 2 . When the boundary condition is of the 
form (8), the condition is trivially satisfied for any t as 
&3 = d = in this case. For the boundary condition 
(9), we have A = bi/b 3 = —b 2 /bs, and the long time 
approximation holds if 



2A 



1^72 



< t. 



(13) 



For the membranes used in our experiments A is of order 
10~ 2 [mm/s] and assuming that we deal with the normal 
diffusion D is roughly 10~ 5 [mm 2 /s]. Thus, the l.h.s. of 
Eq. (13) is estimated as 2 s. Since 10 s is the time step of 
our measurements which extend to 2500 s, the condition 
(13) is fulfilled. We have also checked the condition (13) 
a posteriori, using the values of a and D a obtained by 
means of our method. The l.h.s. of Eq. (13) is again 
about 2 s. 

Let us now discuss the temporal evolution of near- 
membrane layers in the long time approximation. Substi- 
tuting the solution (12) into Eq. (2), we get the equation 
which simplifies to 



H 



;?(( 



1 1 
^ 



One observes that Eq. (14) is solved by 
5{t) = A(a,D a ,K) t a/2 . 



net 
~2~ 



(14) 



(15) 



The near-membrane layer (15) does not depend on the 
parameters b\ and b 2 while the coefficient A can be re- 
calculated into the diffusion constant D a as 



D„ = 



A 2 



2 ) 



(16) 



We have also studied the near-membrane layers be- 
yond the long time approximation using the boundary 
conditions (8) and (9). The condition (8) allows for the 
analytic treatment of Eq. (3) and the solution is of the 
form (12) with 2b\/{b\ — b 2 ) replaced by 1 — a. Thus, 
the formulas derived in the long time approximation are 
exact for Eq (8). When Eq. (9) is used as the boundary 
condition, the solution of subdiffusion equation (3) for 
x > is 
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C(x,t) 



OO 
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71=0 



Co 
a 



(17) 
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;? ( 
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The solutions (12) and (17), which for normal diffusion 
have been discussed in [10,11], qualitatively differ from 
each other but the differences are evident only for times 
which are significantly longer than those studied here. 

Since the formula (17) is analytically intractable we 
have found the time evolution of near-membrane layer 
numerically. As discussed in detail in [5], we have not 
found any difference between the near-membrane layer 
obtained for the concentration profile with the boundary 
condition (8) and with the boundary condition (9). 

Fitting the experimental S(t) by the function At 1 , we 
have found the index a = 2j = 0.90 ± 0.01. It docs 
not much differ from unity but it signals subdiffusion 
due to the small error [5]. With the numerical values 
of inverse Fox functions, we recalculate the coefficient A 
into D a by means of the relation (16). Thus, we get 
D .90 = (9.8 ± 1.0) x 10~ 4 [mm 2 /s - 90 ] for glucose and 
D .go = (6.3 ±0.9) x 10~ 4 [mm 2 /s - 90 ] for sucrose. 
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FIG. 3. The experimentally measured 8 divided by the co- 
efficient A from Eq. (15). The symbols are assigned as in 
Fig. 2 and the line represents the function t 45 . For clarity 
of the plot the error bars are not shown. 

To be sure that Eq. (15), which is used to evaluate D a , 
properly describes the experimental 6(t), we have checked 
the scaling of S(t) suggested by Eq. (15). In Fig. 3 we 
plot the rescaled near-membrane layer 6'(t) — S(t)/A, 
with A from Eq. (15), for all values of k, for glucose and 
for sucrose. The experimental points are represented as 
in Fig. 2. As seen, our experimental data are very well 
described by the function t 0Aa . 

Our method to determine the parameters of subdiffu- 
sion relics on the near-membrane layers. One may ask 
why a and D a are not extracted directly form the con- 



centration profiles which are measured. There are three 
reasons to choose the near-membrane layers: experimen- 
tal, theoretical and practical: 1) Measurement of 5 does 
not suffer from the sizable (~ 10-15%) systematic er- 
ror of absolute normalization of C, as only the relative 
concentration matters for 5. 2) Computed concentration 
profiles depend on the adopted boundary condition at 
a membrane while the condition is not well established 
even for the normal diffusion. The near-membrane layer 
appears to be free of this dependence. 3) When C is fit- 
ted by a solution of the subdiffusion equation, there are 
three free parameters: a, D a and the parameter char- 
acterizing the membrane permeability. Because these fit 
parameters are correlated with each other it is very dif- 
ficult to get their unique values. When S is studied the 
membrane parameter drops out entirely, a is controlled 
by the time dependence of 5(t) while D a is provided by 
the coefficient A. 

The membrane plays only an auxiliary role in our 
method to measure the subdiffusion parameters but the 
transport in membrane systems is of interest in several 
fields of technology [12], where the membranes are used 
as filters, and biophysics [13], where the membrane trans- 
port plays a crucial role in the cell physiology. The dif- 
fusion in a membrane system is also interesting by itself 
as a nontrivial stochastic problem, see e.g. [10]. Thus, 
our study of the subdiffusion in a membrane system, 
which to our best knowledge has not been investigated by 
other authors, opens up a new field of interdisciplinary 
research. 
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